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Abstract 

The  characteristic  parameters  of  the  gas  diffusion-layer  (GDL)  on  cell  performance  and  mass  transfer  of  a  proton  exchange  membrane  fuel 
cell  have  been  investigated  numerically.  A  two-dimensional,  isothermal  and  multi-phase  numerical  model  has  been  established  to  investigate  the 
influence  of  the  GDL  parameters  on  the  transport  phenomenon  and  cell  performance  of  PEM  fuel  cells.  The  porosity  and  thickness  of  the  GDL 
are  employed  in  the  analysis  as  the  parameters.  In  addition,  the  effects  of  liquid  water  and  the  flow  direction  of  the  fuel  and  air  on  the  performance 
are  also  considered  in  this  paper.  The  results  show  that  both  the  porosity  and  thickness  of  the  GDL  affect  the  fuel  cell  performance  significantly, 
especially  the  water  mass  transfer.  It  is  shown  that  the  cell  performance  with  consideration  of  a  liquid  water  effect  is  always  less  than  that  without 
consideration  of  the  liquid  water  effect.  In  addition,  the  cell  performance  with  a  co-flow  pattern  of  fuel  and  air  is  better  than  that  with  a  counter 
flow  pattern. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  is  a 
promising  alternative  power  plant  for  transportation  due  high 
efficiency,  low  emission,  low  noise  and  a  low  operating  tem¬ 
perature  [1-4].  The  gas  diffusion  layers  (GDLs)  of  the  fuel  cell 
is  required  to  provide  both  reactant  gases  to  the  catalyst  layer 
and  removal  of  water  in  either  vapor  or  liquid  form  in  a  typi¬ 
cal  PEMFC.  Although  it  seems  to  be  a  minor  component  in  a 
fuel  cell,  the  GDL  is  one  of  the  most  important  parts  of  a  PEM 
fuel  cell.  A  detailed  study  was  done  by  Moreira  et  al.  [5]  on  the 
influence  of  the  hydrophobic  material  content  in  the  gas  diffusion 
electrode  on  the  performance  of  the  membrane  electrode  assem¬ 
bly  (ME A).  Jordan  et  al.  [6,7]  experimentally  examined  the  influ¬ 
ence  of  the  diffusion-layer  morphology  on  cell  performance. 
They  developed  a  model  of  the  hydrophobilicity  and  porosity  of 
the  diffusion  layer  to  explain  the  influence  of  the  diffusion-layer 
morphology  and  showed  that  a  low-porosity  acetylene  black 
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enhances  water  removal  from  the  ME  A.  Recently,  the  effects 
of  the  fabrication  method  and  the  thickness  of  the  GDL  on  the 
cell  performance  were  experimentally  studied  by  Lee  et  al.  [8]. 
However,  their  study  was  focused  on  the  fabrication  method  of 
the  GDL — rolling,  spraying  and  screen  printing,  in  affecting  the 
cell  performance. 

Because  experimental  work  is  costly,  numerical  modeling 
becomes  an  efficient  and  convenient  approach  to  fuel  cell  anal¬ 
ysis.  For  the  last  decade,  much  effort  has  been  involved  in 
the  development  of  a  numerical  model  with  increasingly  less 
restrictive  assumptions  and  more  physical  complexities.  Several 
examples  of  analysis  of  PEM  fuel  cells  can  be  found  in  the  litera¬ 
ture  [9-12].  The  first  ID  model  of  a  PEM  fuel  cell  was  developed 
by  Springer  et  al.  [13].  Bernardi  and  Verbrugge  [14]  developed 
a  ID  hydraulic  model  and  assumed  that  the  membrane  is  fully 
saturated  with  water  and  that  most  of  the  water  is  transported 
through  the  electrodes  in  the  liquid  phase.  The  first  quasi-2D, 
along-the-channel  model  of  a  PEM  fuel  cell  was  established  by 
Fuller  and  Newmann  [15]  with  the  assumption  of  constant  dif- 
fusivity  of  water  in  the  membrane.  Gurau  et  al.  [16]  presented 
a  comprehensive  model  for  the  entire  sandwich  of  a  PEMFC 
including  the  gas  channels  and  considered  the  gas-liquid  phases 
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Nomenclature 

a  chemical  activity  of  water  vapor  in  cathode 

A  jo  exchange  current  density  (A  m-3) 

C  concentration 

Cp  quadratic  drag  factor 

D  diffusivity  (m2  s  - 1 ) 

F  Faraday  constant,  96485  c  mol-1 

i  current  density  (A  m-2) 

j  current  density  (A  m-3) 

k  permeability  (m2) 

M  molecular  weight 

P  pressure  (atm) 

Pi  partial  pressure  for  i  species  (atm) 

R  universal  gas  constant  (8.314  J  mol-1  K-1) 
s  saturation,  the  ratio  of  the  volume  of  pore  occu¬ 

pied  by  liquid  water  to  the  volume  of  pore  in  the 
porous  medium 

S  source  term  in  momentum  equation  (m  s-2) 

Sc  source  term  of  chemical  reaction  in  the  species 
concentration  equation  (s-1) 

Sj  source  term  in  phase  potential  equation  (A  m-3) 
Sl  source  term  with  consideration  of  liquid  water  in 
the  species  concentration  equation  (s-1) 
t  thickness  (m) 

T  temperature  (K) 

U,  V  velocities  in  the  X  and  Y direction  (ms-1) 

X,  Y  rectangular  coordinate  system  (m) 

Z  number  of  electrons  transferred 

Zf  charge  transfer  coefficient 

Greek  symbols 
a  charge  transfer  rate 

s  porosity 

0  overpotential 

v  kinematic  viscosity  (m2  s-1) 

p  density  (kg  m- 3 ) 

o  electric  conductivity  (1  m-1) 

r  tortuosity  of  the  pore  in  the  porous  medium 

0  membrane  potential 

Superscripts  and  subscripts 

a  quantity  in  anode 

c  quantity  in  cathode 

eff  effective  value 

g  of  gas  diffuser  layer 

H+  for  proton 

H2  for  hydrogen 

H2O  for  water 

i  for  i  species 

m  of  membrane 

O2  for  oxygen 

ref  reference  value 

sat  saturation  pressure  for  water  vapor 

v  in  the  X-direction 

y  in  the  F-direction 


in  the  separate  computation  domains  for  transport  in  the  gas 
distribution  channels.  Hsing  and  Futerko  [17]  developed  a  2D 
model  of  coupled  fluid  flow,  mass  transport  and  electrochem¬ 
istry  of  a  PEMFC  by  taking  into  account  the  dependence  of  the 
diffusion  coefficient  of  liquid  water  in  the  membrane.  However, 
these  2D  models  do  not  resolve  the  catalyst  layers  and  hence 
ignore  the  influence  of  spatial  non-uniformity  of  water  content 
on  the  catalyst  layer  performance. 

Gurau  et  al.  [18]  proposed  a  ID  mathematical  model  to  obtain 
an  analytic  solution  of  the  mass  transport  of  reactant  gas  in  a  half¬ 
cell,  in  which  the  effects  of  the  porosity  and  the  tortuosity  of  the 
GDL  and  the  catalyst  layer  were  explored  due  to  the  fact  that  the 
pores  may  be  partially  filled  with  liquid  water.  Chu  et  al.  [19]  also 
used  a  half-cell  model  to  investigate  the  effects  of  non-uniform 
porosity  on  fuel  cell  performance  in  terms  of  physical  parameters 
such  as  oxygen  consumption,  current  density  and  power  density, 
etc.  Later,  Yan  et  al.  [20]  conducted  a  numerical  study  with  2D 
half-cell  mass  transport  model  to  study  the  effects  of  fuel  channel 
width  and  GDL  porosity  on  the  cross-cell  transport  of  reactant 
gas  and  the  performance  of  a  PEM  fuel  cell.  Recently,  Weber 
and  Newman  [21]  reviewed  on  the  modeling  transport  in  PEM 
fuel  cells  to  discuss  the  different  fuel-cell  models  with  various 
types  of  transport  in  fuel  cells  focusing  on  the  transport  of  the 
various  species  within  the  fuel  cell.  In  this  paper,  they  had  made 
an  overview  for  the  numerical  modeling  work. 

To  consider  the  water  management  of  the  fuel  cell,  many 
of  the  studies  investigated  the  water  transport  of  the  two-phase 
flow  system  in  which  relevant  heat  transfer  problems  were  also 
taken  into  account.  Some  studies  considered  the  physical  domain 
including  the  gas  diffusion  layer,  the  catalyst  layer  and  the  mem¬ 
brane  [22-28],  while  some  focused  on  the  membrane  [29,30].  A 
common  conclusion  from  these  studies  was  that  the  operational 
current  density,  the  humidification  parameter  and  the  membrane 
thickness  had  a  significant  influence  on  the  water  transportation 
in  the  fuel  cell.  From  the  literature  review  above,  it  is  found 
that  a  full-cell  mathematical  model  with  consideration  of  the 
catalyst  layer  and  membrane  water  transport  in  a  PEMFC  has 
not  been  well  examined  yet.  This  motivates  the  present  study. 
The  objective  of  this  work  is  to  establish  a  2D,  full-cell  math¬ 
ematical  model  with  consideration  of  water  saturation,  in  order 
to  investigate  the  effects  of  both  the  porosity  and  size  scale  of 
gas  diffuser-layer  on  the  cell  performance  of  PEMFCs.  Addi¬ 
tionally,  the  effect  of  flow  directions  of  fuel  and  air  on  the  cell 
performance  is  also  taken  into  account  in  the  analysis. 

2.  Analysis 

The  PEM  fuel  cell  model  described  in  this  study  is  developed 
to  analyze  the  characteristics  of  the  GDL.  The  computational 
domain  is  the  full  cell,  which  includes  a  membrane  sandwiched 
between  two  gas  diffusion  electrodes,  and  the  flow  channels  of 
both  the  anode  and  cathode.  Fig.  1  shows  a  schematic  illustration 
of  a  PEMFC  in  the  co-flow  mode  with  the  coordinate  system.  The 
U  and  V  are  the  velocity  components  in  the  X-  and  F-directions, 
respectively.  The  PEMFC  in  the  counter-flow  mode  is  also  exam¬ 
ined  to  study  the  performance  of  a  PEMFC  between  the  co-flow 
and  the  counter-flow  modes  for  comparison.  Consequently,  to 
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Fig.  1.  Schematic  diagram  of  a  complete  PEMFC  and  coordinate  system. 


simplify  the  problem,  a  steady  state,  2D,  multi-species,  and 
along-the-channel  model  of  a  full-cell  PEMFC  is  employed  for 
the  study.  There  are  four  species:  hydrogen,  oxygen,  nitrogen, 
and  water  vapor  considered  in  this  analysis.  Stationary  condi¬ 
tions  are  assumed  in  this  fuel  cell,  also  the  effect  of  gravity  is 
neglected.  It  is  assumed  that  the  electrochemical  reactions  takes 
place  only  in  the  catalyst  layer,  and  the  gas  mixtures  in  the  flow 
channels  are  also  considered  to  be  perfect  gases.  Based  on  the 
definition  of  the  Reynolds  number  and  the  velocity  used  in  this 
work,  the  flow  in  the  fuel  cell  is  laminar.  Therefore,  all  the  trans¬ 
port  equations  were  formulated  for  laminar  behavior.  The  GDL, 
catalyst  layer  and  PEM  are  assumed  to  be  isotropically  porous 
materials. 

According  to  the  descriptions  and  assumption  above,  the 
basic  transport  equations  for  the  2D  PEM  fuel  cell  are  given 
as  the  following: 


•  Continuity  equation: 


dU  dV 

dx  +  3y 


•  Momentum  equation: 


dU  dU  1  dP  ( d2U  d2U\ 

dX  dY  pdX  V  9X-2  dY2  J 

dV  dV  1  3 P  ( d2V  d2V\ 

dx  dY  pdY  V dx  dY  J 


+  Sx 


H"  Sy 


•  Species  equation: 


£&,eff 


+Sc  +  SL 


(2) 

(3) 


In  the  momentum  equations,  Sx  and  Sy  stand  for  the  source 
terms  based  on  the  Darcy’s  drag  forces  in  the  X  and  Y  directions 
imposed  by  the  pore  walls  on  the  fluid,  and  usually  cause  in  a 
significant  pressure  drop  across  the  porous  media.  The  details 
of  Sx  and  Sy  for  different  layers  are  listed  in  Table  1 .  In  Table  1 , 
the  £*:,e ff  is  the  effective  porosity,  r  the  tortuosity  of  the  pores 
in  porous  medium,  Cp  represents  the  quadratic  drag  factor,  k 
the  permeability  for  the  porous  medium,  Zf  the  charge  transfer 
coefficient,  CH+  is  the  concentration  of  proton,  F  is  the  Faraday 
constant  and  0  represents  the  membrane  potential.  In  the  anal¬ 
ysis,  Blake-Kozeny  equation  [31]  is  used  to  model  k  as  below: 


|A2p\ 

~  c3  1 

V150/ 

(1  -  ££,e ff)2 

where  Dip  =  6Rys,  and  is  the  volume- to- surface  ratio  of  the 

porous  material.  The  parameters  ja  and  jc  in  Table  1  indicate  the 
current  density  at  the  anode  and  cathode  sides,  respectively,  and 
can  be  described  by  the  following  Butler- Volmer  equations: 


CH2 

z^ref 

CH2 


1/2  r 


t(aaF/RT)r] 


1 

Q(acF/RT)r] 


e(aa  F/RT)r) 


1 

Q(acF/RT)ri 


where  A  is  the  exchange  current  density,  aa  and  ac  the  electric 
charge  transport  rates  in  anode  and  cathode  catalyst  layers,  rj  the 
over-potential,  R  the  gas  constant  and  T  is  the  temperature  of  the 
fuel  cell. 


Table  1 

Detailed  expressions  of  the  source  terms  in  the  governing  equations 


Gas  diffuser  layer  —  V£^eff  U - g,ef^- —  s/u1  +  V2 

2 

Catalyst  layer  -  U  -  VC/2  +  V2 

vs2  ,,  £3  ffCFPU  , - = - ^ 

PEM  - jH.rff  - Vf/2  +  V2  + 

;Zf  Ch+F-V0-§ 


y_  ei-‘fPV  ^tj2  +  y2 
k  Vk 

-  t;£c.eff  y  _  Jjj2  y2 

vs2  ....  e3  ffCppV  , — ^ ^ 

_  m.effy  _  m,eff  ^  ^  jj2  +  y2  + 

k-ZfCn+F.V0.^ 


Di,e  ff 

T 

0 

Di  '  £g,eff 

1.5 

H2:  —  (l/2FCa)ja; 
H20:  (1/2 FCc)jc 

02: 

-(1/4  FCc)jc- 

Di  '  £c,eff 

1.5 

ff  A',eff,H+ Gt+  ( 

d2& 

dx2 

+  d2(p ^ 

+  dY 2  J 

Di  '  £m,eff 

6 
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In  Eq.  (4),  Sc  is  the  production  rates  of  i- th  species  in  gas 
phase,  and  A  ,eff  is  the  effective  mass  diffusivity,  and  Sfi  repre¬ 
sents  the  quality  of  liquid  water  in  order  to  investigate  the  liquid 
water  effect  in  this  work.  The  model  modifies  the  mass  diffusiv¬ 
ity  due  to  the  liquid  water  filling  the  pores  in  the  porous  media 
and  the  liquid  generation  in  the  species  equation.  When  the  par¬ 
tial  pressure  of  water  vapor  is  greater  than  the  saturation  pressure 
of  water  vapor,  the  water  vapor  is  assumed  to  condense  and  fill 
the  pore  in  the  porous  media.  Therefore,  Sl  can  be  evaluated  by 
the  following  [6] : 


^H20^c(£eff^H2o/ PRT)(Pr20 
ke£effS(Ps^  -  Ph2o), 


^sat)>  if  ft[20  >  ^sat 

if  ^H20  <  ^sat 

(8) 


where  the  M  is  the  molecular  weight  and  kc  and  ke  are  the 
condensation  and  evaporation  rate  constants,  respectively.  The 
saturation  pressure  of  water  can  be  expressed  as  [13] 

psat  =  i0_2-1794+0-02953r_9-1837xl0_5:r2+L4454xl0_7:r3  (9) 


In  addition,  the  saturation,  s,  is  defined  as  the  ratio  of  the 
volume  of  pore  occupied  by  liquid  water  to  the  volume  of  pore  in 
the  porous  medium,  then  the  effective  porosity  of  porous  media 
is  modified  to  account  the  liquid  water  effect, 

£eff  =  e(l  -  s )  (10) 


In  order  to  evaluate  the  distributions  of  the  local  current  den¬ 
sity,  the  phase  potential  equation  should  be  solved, 


9  f  d&\  d  f  _ 

3x  Vm9X )  +  dY  Vm  dY  )  -  Sj 


(11) 


dix  diy 

- 1 - -  =  jc  at  cathode  (18) 

dX  dY  J 

Boundary  conditions  for  the  dependent  variables  of  the  trans¬ 
port  equations  at  the  interfaces  between  different  layers  of  the 
same  domain  are  not  required.  A  fully  developed  flow  condition 
is  assumed  at  the  channel  outlet,  therefore, 


(19) 


The  boundary  conditions  at  the  gas  flow  channel  walls  are 


dCi 

U  =  V  =  — 1=0  (20) 

dX 

In  practical  situations,  the  physical  properties,  such  as  veloc¬ 
ity  and  concentration,  and  their  gradients  are  continuous  on  the 
interface.  So  the  natural  boundary  conditions  on  the  interface  are 
the  same  velocity,  same  concentration  and  the  same  gradients  of 
both  velocity  and  concentration.  At  the  interfaces  between  the 
gas  diffusers  and  the  gas  channels,  the  following  boundary  con¬ 
ditions  are  used, 


dv 

£eff,x+  dx 

x=x+ 

dV 

dX 

5 

x=x- 

II 

X 

+ 

II 

II 

X 

1 

(21) 

dCi 

£eff,x+  dx 

dCi 

dX 

Q,X=X+  —  Q,X=X- 

(22) 

Similar  conditions  are  employed  for  the  interfaces  between 
the  gas  diffuser  layer  and  the  catalyst  layers  and  the  interfaces 
between  the  catalyst  layers  and  membrane  can  be  expressed  as 
follows: 


where  Sj  is  — ya  in  the  anode,  —jc  in  the  cathode  and  0  in  the  mem¬ 
brane;  crm  is  the  electric  conductivity  of  the  membrane  which 
can  be  calculated  by  the  equation  developed  by  Springer  et  al. 
[13]: 


£&ff,X+ 


dv 

dx 


dv 

-  £eff,X~  ™ 

x=x+  dX 


Vx=x+  —  Vx=x~ 


X=X- 


(23) 


crm(T)  =  ovf  exp 


1268 


1 


303 


1 

T 


(12) 


and  the  reference  electric  conductivity  is 
=  0.005 139A-  0.00326 


(13) 


0.043  +  17.81a  —  39.85a2  +  36.0a3  0  <  a  <  1 

A  =  <  ~  ~  (14) 

14  +  1.4(a  —  1)  1  <  a  <  3 

where  a  is  the  activity  of  water  vapor  at  the  cathode  side.  Using 
the  following  relations  between  the  phase  potential  0  and  current 


density  i: 

d0 

(15) 

dX 

d0 

ly-  qy 

(16) 

Therefore,  Eq.  (11)  can  then  be  reduced  to  be 

dix  dl 


H - 2  —  ja  at  anode 


dX  dY 


(17) 


£&ff,X+ 


dQ 

dX 


dCi 

—  £eff,X~ 

x=x+  dX 


X=X- 


Q,X=X+  —  Ci,X=X~ 

(24) 


The  boundary  conditions  for  the  phase  potential  at  the  inter¬ 
face  between  the  catalyst  layer  and  the  membrane  are  0  =  0  at 
the  anode  side,  and  d0/dX=O  at  the  cathode  side.  Because  the 
phase  potential  is  a  linear  distribution  in  the  membrane,  the  phase 
potential  boundary  condition  can  be  written  as  d0/dY=  0. 


3.  Numerical  method 

The  solution  to  the  governing  equations  is  performed  using 
a  finite  volume  scheme  by  dividing  the  model  domain  into  a 
number  of  cells  as  control  volumes.  In  the  finite  volume  method, 
the  governing  equations  are  numerically  integrated  over  each  of 
these  computational  cells  or  control  volumes.  The  finite  volume 
method  exploits  a  collocated  cell-centered  variable  arrangement, 
that  implies  all  the  dependent  variables  and  material  properties 
are  stored  at  the  cell  center.  The  average  value  of  any  quantity 
within  a  control  volume  is  given  by  its  value  at  the  cell  center. 
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The  governing  equations  can  be  expressed  in  the  form  of  a 
generalized  transport  equation: 

V  •  (pucp  -  L0V0)  =  (25) 

where  0  denotes  the  general  dependent  variable,  r $  the 
exchange  coefficient,  S $  the  source  term,  u  represents  the  veloc¬ 
ity  vector,  and  p  is  the  density.  With  the  discretization  of  the 
governing  equations,  the  coupled  finite-difference  equations  can 
be  expressed  in  the  form  of 

<2p0p  =  <7E0E  +  <7W0W  +  <7N0N  +  <2s0S  +  Sf  (26) 

where  0p  is  the  value  of  0  at  the  current  point  P,  0e-  .  .0s  stand 
for  the  values  of  the  grid  points  adjacent  to  the  point  P,  and 
cip. .  .as  are  known  as  the  link  coefficients. 

In  this  work,  the  non-uniform  grid  system  of  85  x  70  is 
employed  for  the  analysis.  In  order  to  examine  the  grid  inde¬ 
pendence  of  the  predictions,  coarse  and  fine  grid  systems  are 
considered  in  the  preliminary  tests.  Effects  of  the  grid  number 
on  the  predictions  of  local  current  density  are  shown  in  Fig.  2. 
The  maximum  deviations  among  the  computations  on  the  grids 
of  50  x  35,  85  x  70  and  130  x  140  are  less  than  3%.  Therefore, 
the  grid  system  of  85  x  70  points  seems  to  be  sufficient  to  resolve 
the  behaviors  of  the  reactant  gas  transport  in  the  present  PEMFC 
model. 

4.  Results  and  discussions 

The  inlet  conditions  for  the  PEM  fuel  cell  are  the  inlet  pres¬ 
sure  with  101.3  kPa,  inlet  temperature  with  333.15  K,  relative 
humidity  with  100%  and  inlet  velocity  with  1ms-1  for  both 
anode  and  cathode.  The  total  length  of  the  flow  channel  is  14  cm, 
and  the  cross-section  of  the  flow  channel  is  1  mm  x  1  mm.  The 
thicknesses  of  the  catalyst  layer  and  the  membrane  are  fixed 
and  taken  to  be  0.0000287  and  0.00175  m,  respectively.  In  this 
study,  the  effects  of  the  GDL  porosity  ranging  from  0.3  to  0.6 
and  thickness  ranging  from  0.0002  to  0.0006  m  are  examined. 
Results  without  liquid  water  are  also  investigated  as  well  as  the 
flow  direction  effect. 


Fig.  3.  Effects  of  the  GDL  porosity  on  the  polarization  curves  of  the  PEMFC 
with/without  liquid  water  effects. 

Fig.  3  shows  the  polarization  ( I-V)  curves  with  various  GDL 
porosity  of  a  fuel  cell  to  investigate  the  influence  of  GDL  poros¬ 
ity  on  the  cell  performance.  The  results  without  liquid  water 
effects  are  also  presented  for  comparison.  It  shows  in  Fig.  3  that 
the  effect  of  the  GDL  porosity  on  the  cell  performance  is  sig¬ 
nificant  at  low  operating  conditions.  However,  at  high  operating 
conditions,  the  influence  of  GDL  porosity  on  the  I-V  curve  is 
negligibly  small.  It  is  also  observed  that  the  cell  performance  is 
increased  as  the  GDL  porosity  increases  with/without  consid¬ 
eration  of  liquid  water  effect.  For  a  GDL  with  higher  porosity 
contains  much  void,  which  allows  more  gas  reactant  transfer 
into  catalyst  layer,  in  turns,  more  chemical  reaction  occurs  and 
resulted  in  a  better  cell  performance.  In  addition,  it  is  noticed 
that  the  I-V  curve  is  over-predicted  when  the  liquid  water  effect 
is  not  taken  into  account  in  the  modeling.  This  can  be  explained 
by  the  fact  that  the  void  in  the  GDL  is  filled  with  the  liquid  water, 
which  in  turn,  causes  the  reduction  of  mass  transfer.  It  is  also 
found  in  Fig.  3  that  the  cell  performance  differences  between  the 
results  with  and  without  consideration  of  liquid  water  effects  are 
small  at  high  operating  voltage  conditions.  It  means  that  the  fuel 
transport  in  the  PEMFC  can  be  treated  as  single  gas  phase  at 
high  voltage  conditions.  However,  the  liquid  water  effects  on 
the  cell  performance  are  remarkable  and  cannot  be  neglected 
in  the  modeling  at  low  voltage  conditions.  This  confirms  the 
fact  that  the  mass  transports  are  significant  at  lower  voltage 
operating  conditions  and,  in  turn,  more  water  is  generated  in 
the  catalyst  layer  of  the  cathode  side.  Therefore,  two-phase 
flow  effects  should  be  considered  under  low  operating  voltage 
conditions. 

It  is  important  in  the  design  of  a  PEM  fuel  cell  to  realize  the 
distribution  of  the  fuel  gases  in  the  catalyst  layer.  The  hydro¬ 
gen  mass  flux  distributions  in  the  catalyst  layer  at  the  operating 
voltage  of  0.2  V  with  liquid  water  effect  are  shown  in  Fig.  4 
for  the  GDL  porosity  is  0.3,  0.4,  0.5  and  0.6.  The  interface  of 
catalyst  layer  and  membrane  is  on  the  right  side,  while  the  inter¬ 
face  of  the  GDL  and  the  catalyst  layer  is  on  the  left.  It  is  seen 
that  the  mass  flux  of  hydrogen  decreases  along  the  axial  direc¬ 
tion  due  to  the  reaction  of  the  hydrogen.  Besides,  the  hydrogen 
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Fig.  4.  Effects  of  the  GDL  porosity  on  the  fuel  mass  flux  in  the  anode  catalyst  layer:  (a)  £g  =  0.3;  (b)  £g  =  0.4;  (c)  £g  =  0.5;  (d)  £g  =  0.6. 


Saturation 


Fig.  5.  Effect  of  the  GDL  porosity  on  the  liquid  water  distribution  for  a  full-cell  PEMFC:  (a)  £g  =  0.3;  (b)  £g  =  0.4;  (c)  £g  =  0.5;  (d)  £g  =  0.6. 
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mass  flux  decreases  along  the  X-direction.  Results  show  that  the 
mass  flux  of  the  hydrogen  decreases  as  the  porosity  increases. 
This  means  that  the  consumption  of  the  fuel  gases  increases 
as  the  GDL  porosity  increases.  This  again  confirms  the  result 
mentioned  above.  The  liquid  water  saturation  distribution  at  the 
operating  voltage  of  0.2  V  for  the  whole  fuel  cell  is  shown  in 
Fig.  5  with  different  GDL  porosities.  As  saturation  equals  to 
0  indicates  that  there  is  no  liquid  water  produced,  while  sat¬ 
uration  equals  to  1  indicates  that  the  porous  material  is  filled 
fully  with  liquid  water.  An  overall  inspection  on  Fig.  5  indicates 
that  the  GDL  porosity  affects  on  the  liquid  water  production 
dramatically.  This  may  be  explained  by  the  fact  that  the  electro¬ 
chemical  reaction  rate  is  high  at  lower  operating  voltages  and 
more  mass  transfer  for  higher  GDL  porosity,  which  can  thus 
generate  more  liquid  water.  It  is  also  seen  that  the  liquid  water 
increases  along  the  axial  location  and  the  peak  value  of  liquid 
water  locates  at  the  cathode  catalyst  layer.  This  is  due  to  the  liq¬ 
uid  water  generation  occurring  at  the  catalyst  layer  of  the  cathode 
side. 

The  effect  of  the  GDL  thickness  on  the  I-V  curve  of  a  PEM 
fuel  cell  is  shown  in  Fig.  6.  It  is  obvious  that  the  cell  performance 


Fig.  6.  Effects  of  the  GDL  thickness  on  the  polarization  curves  of  the  PEMFC 
with/without  liquid  water  effects. 
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Fig.  7.  Effects  of  the  GDL  thickness  on  the  fuel  mass  flux  in  the  anode  catalyst  layer:  (a)  tg  =  0.0002  m;  (b)  tg  =  0.0004  m;  (c)  tg  =  0.0006  m. 
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with  or  without  liquid  water  effect  increases  as  the  GDL  thick¬ 
ness  decreases,  especially  at  lower  operating  voltage  conditions. 
For  a  thinner  GDL,  the  concentration  gradient  becomes  larger 
causing  a  higher  mass  transfer  into  the  catalyst  layer.  Therefore,  a 
larger  current  density  takes  place  at  a  lower  operating  condition. 
However,  it  does  not  indicate  the  same  trend  at  higher  operat¬ 
ing  conditions  for  V>0.6  V.  Fig.  7  presents  the  hydrogen  mass 
flux  distribution  in  the  catalyst  layer  with  various  GDL  thick¬ 
nesses  at  the  operating  voltage  of  0.2  V.  It  is  observed  that  the 
hydrogen  mass  flux  increases  as  the  GDL  thickness  increases. 

Saturation 

0  0.5  • 


(a) 


(b) 


(c) 

Fig.  8.  Effect  of  the  GDL  thickness  on  the  liquid  water  distribution  for  a  full-cell 
PEMFC:  (a)  tg  =  0.0002  m;  (b)  tg  =  0.0004  m;  (c)  tg  =  0.0006  m. 


Fig.  9.  Effects  of  the  flow  direction  of  the  reactant  gases  on  the  polarization 
curves  of  the  PEMFC. 

This  is  due  to  the  fact  that  the  concentration  gradient  is  smaller 
with  a  thicker  GDL,  which  results  in  a  smaller  mass  transfer  of 
fuel  gases.  The  liquid  water  saturation  distributions  with  various 
GDL  thicknesses  for  the  whole  cell  at  the  operating  voltage  of 
0.2  V  are  presented  in  Fig.  8.  The  saturation  increases  with  a 
decrease  in  the  GDL  thickness.  This  is  because  more  fuel  gas 
transfer  into  the  catalyst  layer,  which  results  in  more  liquid  water 
generation. 

The  effect  of  the  flow  direction  of  reactant  fuels  in  both  anode 
and  cathode  on  the  cell  performance  is  also  studied  in  this  work. 
Fig.  9  presents  the  effects  of  the  GDL  porosity  on  the  cell  per¬ 
formance  with  the  reactant  fuels  in  anode  and  cathode  being 
co-flow  and  counter-flow  conditions.  It  is  found  that  the  cell 
performance  increases  as  the  GDL  porosity  increases  for  both 
co-flow  and  counter-flow  conditions.  It  is  also  discovered  that 
the  cell  performance  with  counter-flow  condition  is  better  at 
higher  operating  voltage  conditions,  while  the  cell  performance 
with  co-flow  condition  is  better  at  lower  operating  voltage  con¬ 
ditions.  Fig.  10  depicts  the  distributions  of  hydrogen  mass  flux 
in  the  catalyst  layer  with  counter-flow  at  the  operating  voltage 
of  0.2  V.  It  is  observed  that  the  hydrogen  mass  flux  decreases 
as  the  GDL  porosity  increases  with  counter-flow  condition.  The 
reason  is  the  same  as  that  with  co-flow  condition.  Comparison 
of  the  corresponding  results  in  Figs.  4  and  10  discloses  that  the 
hydrogen  mass  flux  in  counter-flow  condition  is  lower  than  that 
in  co-flow  condition  for  various  GDL  porosities.  This  indicates 
that  the  hydrogen  consumption  is  smaller  with  counter-flow  con¬ 
dition.  This  is  consistent  to  the  result  that  the  cell  performance 
of  co-flow  condition  is  better  than  that  of  counter-flow  condition 
at  lower  operating  voltages.  The  liquid  water  saturation  dis¬ 
tributions  at  the  operating  voltage  of  0.2  V  with  counter-flow 
condition  in  a  full  cell  for  various  GDL  porosities  are  shown  in 
Fig.  1 1.  It  is  found  that  the  cathode  GDL  and  catalyst  layer  are 
also  filled  with  liquid  water.  Comparison  of  the  corresponding 
curves  in  Figs.  5  and  11  indicates  that  the  more  liquid  water 
is  generated  at  the  end  of  the  cell  with  counter-flow  condition, 
especially  with  £g  =  0.6. 
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Fig.  10.  Effects  of  the  GDL  porosity  on  the  fuel  mass  flux  in  the  anode  catalyst  layer  with  counterflow:  (a)  sg  =  0.3;  (b)  eg  =  0.4;  (c)  eg  =  0.5;  (d)  eg  =  0.6. 


Saturation 


Fig.  11.  Effect  of  the  GDL  porosity  on  the  liquid  water  distribution  for  a  full-cell  PEMFC  with  counterflow:  (a)  sg  =  0.3;  (b)  eg  =  0.4;  (c)  sg  =  0.5;  (d)  eg  =  0.6. 
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5.  Conclusion 

A  full-cell  numerical  model  with  steady,  2D,  isothermal, 
multi- species  has  been  developed  to  investigate  the  influence 
of  the  gas  diffusion-layer  parameters  on  the  PEMFC.  The  liq¬ 
uid  water  effect  is  also  considered  in  this  study.  The  effects  of 
the  GDL  porosity,  the  GDL  thickness  and  the  flow  direction 
of  the  reactant  gas  on  the  mass  transport  and  cell  performance 
are  examined  in  detail.  A  summary  of  the  major  results  are  as 
follows: 

1 .  It  is  found  that  the  mass  transport  for  both  the  fuel  and  air 
increases  as  the  porosity  of  the  GDL  ranging  from  0.3  to  0.6 
increases.  Therefore,  more  reactant  gases  transfer  into  the 
catalyst  layer,  which  in  turn,  leads  to  more  chemical  reaction 
and  more  reactant  gases  consumed.  This  results  in  a  better 
performance  for  the  fuel  cell  with  a  higher  GDL  porosity 
within  the  range  of  0.3-0. 6. 

2.  The  thickness  of  the  GDL  affects  the  cell  significantly, 
especially  at  lower  operating  voltages,  the  cell  performance 
increases  as  the  GDL  thickness  ranging  from  0.0002  to 
0.0006  m  decreases.  This  is  due  to  the  fact  that  higher  con¬ 
centration  gradients  result  from  the  decrease  in  the  GDL 
thickness,  which  in  turn,  results  in  a  higher  mass  transfer 
to  the  catalyst  layer  and  more  electrochemical  reaction. 

3.  Liquid  water  effects  on  the  cell  performance  are  remarkable 
and  cannot  be  neglected  at  lower  operating  conditions.  There¬ 
fore,  two-phase  flow  effects  should  be  considered  at  low 
operating  conditions.  However,  the  cell  performance  differ¬ 
ences  between  the  results  with  and  without  consideration  of 
the  liquid  water  effects  are  small  at  higher  operating  voltage 
conditions. 

4.  The  cell  performance  with  the  counter-flow  condition  is 
larger  at  the  higher  operating  voltage  conditions,  while  the 
cell  performance  with  a  co-flow  condition  is  larger  at  a  lower 
operating  voltage  conditions.  This  is  due  to  the  fact  that  more 
liquid  water  is  generated  under  the  counter-flow  condition  at 
lower  operating  voltages,  and  the  voids  in  the  GDL  are  filled 
with  liquid  water,  which  in  turn,  causes  a  reduction  of  mass 
transfer  and  a  lower  cell  performance. 
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